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ON  BOUNDARY  EXTRAPOLATION  AND  DISSIPATIVE  SCHEMES  FOR  HYPERBOLIC  PROBLEMS 

Moshe  Goldberg* 

Department  of  Mathematics 
University  of  California 
Los  Angeles,  California,  90024 

ABSTRACT.  In  this  note  we  consider  dissipative,  stable  approximations 
to  well-posed  linear  hyperbolic  initial  value  problems  in  the  quarter  plane 
x > 0,  t > 0.  We  show  that  if  boundary  values  are  determined  by  extra- 
polation, then  stability  is  maintained.  This  result  was  first  suggested 
by  Kreiss,  and  proved  explicitly  by  the  author.  The  proof  is  reviewed 
here  using  a stability  criterion  for  a certain  family  of  boundary  conditions 
due  to  Goldberg  and  Tadmor. 

The  Lax-Wendroff  scheme  and  other  dissipative  approximations  are 
applied  to  a test  problem.  As  expected  from  Gustafsson' s rate-of- 
convergence  theory,  computations  verify  that  if  the  boundary  extra- 
polation and  the  difference  scheme  have  equal  order  of  accuracy,  then 
this  order  is  preserved. 

1.  INTRODUCTION.  Consider  the  conservation  law 

(la)  du/St  + dg(u)/cbc  = 0,  x > 0,  t > 0, 
and  assume  that  the  associated  initial  value  prr’ 

(lb)  u(x,0)  = f(x) 

2 

is  well -posed  in  L (0,«),  so  no  boundary  values  are  required  at  x = 0. 
This  assumption  implies  that  characteristic  lines  do  not  carry  information 
from  the  exterior  of  the  domain  x > 0,  t > 0 inward. 

To  approximate  the  initial  value  problem  (l),  we  introduce  a mesh 
size  Ax  > 0,  At  > 0;  a grid  function  vy(t)  = v(vAx,t),  v = 0,+l,+2,...; 
and  a consistent,  explicit  finite  difference  scheme 

(2)  vv(t  + At)  = S{vv_r(t),...,vv+p(t)};  v = 1,2,3,..., 

r,p  being  fixed  integers. 
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Since  nonlinearity  in  (la)  leads  to  nonlinear  dependence  of 
v (t  + At)  on  the  components  of  v (t), . . . .v,^.  (t),  we  are  unable  to 

v V*r  Vrp 

be  more  specific,  at  this  stage,  about  the  structure  of  the  scheme  S. 

p 

However,  we  assume  that  S is  L -stable,  in  case  it  is  applied  to  the 
pure  initial  value  problem  for  -<*>  < x < <». 

Usually,  r > 0,  so  it  is  impossible  to  approximate  (l)  by  (2) 
without  specifying  boundary  values  at  r grid  points  in  some  left 
neighborhood  of  the  boundary  x = 0.  Thus,  we  admit  boundary  conditions 
of  the  form 


(3) 


V*>  ■ £ W’-  + 


where  the  coefficients  c^  and  s > 1 are  fixed.  That  is,  having  the 
values  vv(t),  v > 1,  computed  by  the  basic  scheme  (2),  we  proceed,  at 
each  time  step,  by  using  (3)  to  determine  v^(t),  u = 0,-1,..., -r  + 1, 
in  that  order. 

A natural  way  to  choose  the  boundary  conditions  in  (3)  would  be  to 
employ  extrapolation  of  degree  s - 1 — a procedure  which  is  of 

accuracy  of  order  s.  More  explicitly,  we  extrapolate  from  vn(t),...,v  (t) 
to  vQ(t) ; then  from  vQ(t), . . .,vs  l(t)  to  v_x(t),  etc.  With  the  use 
of  Stirling's  extrapolation  formula,  (3)  becomes 


(>0 


V‘>  - £ 0)<-« 


3+1 


U — 0,...,— r + 1. 


The  main  purpose  of  this  note  is  to  study  the  influence  of  boundary 
extrapolation  on  the  stability  of  the  numerical  algorithm.  This  question 
is  discussed  in  Section  2,  where  we  consider  a scalar  linear  conservation 
law  which  we  approximate  by  a dissipative  scheme.  In  this  simple  case  it 
is  shown  that  boundary  extrapolation  maintains  stability.  This  result, 
which  was  first  suggested  by  Kreiss,  [5],  was  proven  explicitely  in  [l], 
using  Kreiss'  theory,  [6],  for  dissipative  approximations  of  mixed  initial- 
boundary value  problems. 

In  fact,  the  above  assertion  is  an  immediate  corollary  of  a forth- 
coming work  by  Goldberg  and  Tadmor,  [2],  which  provides  stability  criteria 
for  seme  general  families  of  boundary  conditions,  including  those 
presented  in  (3). 


Finally,  in  Section  3,  the  Lax-Wendroff  scheme,  [7],  and  a new 
5-point  dissipative  approximation  by  Gottlieb  and  Turkel,  [3],  are  applied 
to  a test  problem.  The  numerical  results  support  Gustafsson' s rate-of- 
convergence  theory,  [4],  by  showing  that  the  accuracy  order  of  the  basic 
scheme  is  maintained,  if  the  extrapolation  at  the  boundary  is  of  the 
same  order.  The  computation  were  carried  out  at  the  Campus  Computing 
Network  of  the  University  of  California  at  Los  Angeles. 


2.  STABILITY  ANALYSIS.  From  now  on  we  restrict  attention  to  the 
linear,  scalar  version  of  (l),  namely  to  the  initial  value  problem 

(5)  du/dt  + adu/dx  =0;  a = const.;  x > 0,  t > 0;  u(x,0)  = f(x), 

which  is  well -posed  if  and  only  if  a < 0. 

Our  explicit  approximation  in  (2)  becomes 


(6) 


vy(t  + At)  = Qvv(t),  v = 1,2,3,..., 
« - £_r  % - Yv+1>  r > °> 


where  the  constants  depend  on  a and  on  the  fixed  ratio  X = At/Ax, 

and  initial  values  are  determined  by 

vv(0)  = fv,  v = 1,2,3,...  . 

The  assumption  of  dissipativity  is  that  for  some  & > 0 and  natural 
<*i,  the  amplification  factor  of  the  scheme, 


satisfies 


S(i) 


a eiJ§ 
aje  » 


-7T  < 5 < 7T, 


1 3(0 1 < 1 - &ui2“,  v ui  < tt. 


Thus,  it  is  evident  that  Q(  |)  is  now  power  bounded  (by  l),  which  is 
well  known  to  be  equivalent  to  the  (strong)  stability  of  our  basic  scheme. 
Introducing  the  boundary  conditions  (4),  the  concept  of  stability 


becomes  considerably  more  complicated,  and  we  review  it  briefly.  Let 
H s H(ax)  be  the  space  of  all  grid  functions,  w = {wv)”_  ^ , which 
satisfy  £_r+^|wv|^  < °°  and  fulfill  the  boundary  conditions  in  (4).  If 
inner  product  and  norm  are  defined  by 

CO 

(v,w)  = Ax  s V w , llwl)2  = (w,w), 

V— H-l 

2 

then  H becomes  a discrete  analogue  of  L (0,x). 

Having  constructed  H,  we  realize  that  our  finite  difference 
algorithm  in  (4)  and  (6)  defines  a linear,  bounded  operator,  G : H -* H, 
such  that  the  numerical  solution  v satisfies, 

v(t  + At)  = Gv(t),  for  v(t)  e H. 


Since 


v(t)  = Gmv(0)  for  t = mAt,  m = 1,2,3,..., 

stability  means  that  the  powers  of  G are  uniformly  bounded,  i.e.,  that 
for  some  constant  K, 

||Gm||<K,  m = 1,2,3,...  . 

We  are  now  ready  to  state  the  main  result: 

THEOREM  1.  Let  the  initial  value  problem  (5)  be  approximated  by  an 
arbitrary,  dissipative  (stable)  scheme  of  type  (6),  which  is  complemented 
by  boundary  extrapolation  of  arbitrary  order.  Then,  the  overall  numerical 
algorithm  is  stable. 

The  proof  which  is  laid  out  in  [l],  is  a direct  but  somewhat  lengthy 
application  of  Kreiss1  stability  theory,  [6],  for  dissipative  schemes. 

As  required  by  Kreiss'  criterion,  the  problem  was  to  show  that  the  corres- 
ponding operator  G has  no  eigenvalues  z in  the  unit  disk. 

In  a forthcoming  paper,  [2],  Goldberg  and  Tadmor  use  Kreiss'  theory 
to  provide  a particularly  simple  stability  condition  in  the  case  where 
scheme  (6)  is  augmented  by  boundary  conditions  of  type  (2).  This  condition 
is  rephrased  as  follows: 


THEOREM  2 (Goldberg,  Tadmor).  Let  (6)  be  an  arbitrary,  dissipative 
(stable)  approximation , augmented  by  boundary  conditions  of  type  (2), 
then  the  overall  algorithm  is  stable  if 

E c.K*^  ^ 1,  V*  with  | k | <1. 

J=-r  J 

Theorem  2,  which  is  actually  independent  of  the  basic  scheme,  yields 
Theorem  1 immediately.  For,  considering  the  boundary  conditions  in  (4), 
we  want  to  show  that 

E Qj(-l)^+1  ^ 1 for  k with  |k|  <1, 

i.e.,  that  for  all  k with  |k|  <1, 

(1  - K)s  . E (j)(-«)J  * o. 

The  last  inequality  holds,  and  Theorem  1 follows. 

3.  NUMERICAL  RESULTS.  Consider  the  test  problem 

(7)  du/dt  - du/dx  = 0;  x > 0,  t > 0;  u(x,0)  = sin  2ttx 
whose  analytic  solution  is 

u(x,t)  = sin  2tt(x  + t). 

The  second  order  accurate  Lax-Wendroff  scheme  (L-W),  [7],  is  in  this 

case 

(8)  v(t  + At)  = ^(X  - l)vv-1(t)  + (l  - X2)vv(t)  + |(X  + l^^t),  X = 

and  it  is  well  known  (e.g.,  [8,  Chapter  12])  that  dissipativity,  and 
hence  stability,  are  guaranteed  if  X < 1. 

In  order  to  apply  (8)  to  (7)  we  need  to  specify  only  one  boundary 
value,  vQ(t),  which  according  to  (4),  is  given  by 

T°(t)  ’ Jl(’)('1)3+1 

Here  the  accuracy  of  the  boundary  extrapolation  is  of  orders  s where 
s is  arbitrary. 


In  Table  1 we  compare  the  L-W  results  with  the  analytic  solution 
at  t = 1.  The  H-norm  of  the  error  was  computed  over  the  interval 
0 < x < 1,  and  is  defined  by 


(9) 


J = l/Ax. 


Ax 

s 

m 

116  ”(0,1) 

.05 

2 

40 

5.57  - 2 

.025 

2 

80 

1.39  - 2 

.05 

1 

40 

7.03  - 2 

.025 

1 

80 

2.23  - 2 

Table  1.  L-W  results  at  t = 1; 

X = l/2;  m = t/At  is  number  of  time  steps. 

Gustafsson,  in  his  rate -of -convergence  theory,  [4],  has  discussed 
situations  similar  to  the  one  under  consideration.  He  has  shown  that  in 
order  to  maintain  the  accuracy  of  the  basic  scheme,  it  is  sufficient  to 
employ  boundary  conditions  of  the  same  order  of  accuracy.  Indeed,  Table  1 
suggests  that  L-W' s second  order  accuracy  is  maintained  if  the  boundary 
extrapolation  is  linear  (s  = 2),  but  is  reduced  if  s = 1. 

A second  example  is  concerned  with  a family  of  centered,  5-point, 
dissipative  schemes  by  Gottlieb  and  Turkel  (G-T),  [3].  The  family,  given 
in  (2.4)  of  [3],  depends  on  two  parameters  a and  a.  Choosing  a = l/2, 
= 1,  and  linearizing,  we  obtain  an  approximation  to  (7)  of  the  form 

vv(t  + At)  = - |)vv.2(t)  + X(\  - |)vv_1(t)  + (1  - !>2)vv(t) 

♦ x<x  + §K.iW  - 1(2 ♦ f x ' S> 

where  the  dissipativity  condition  is  X < ^2/2. 

Now  we  need  two  boundary  values  which  are  given  by 

u ■ °<-1- 


The  error-norms  in  Table  2 are  computed,  as  in  the  previous  case. 


Table  2.  G-T  results  for  t = 1. 

Unlike  the  L-W  scheme  which  is  of  second  order  accuracy  both  in 
time  and  space,  the  G-T  approximation  is  of  second  order  in  time  and 
fourth  order  in  space.  Since  the  boundary  extrapolation  is  taken  only 
with  respect  to  the  space  variable,  it  should  be  expected  that  in  order 
maintain  the  fourth-order  accuracy  in  x,  we  have  to  utilize  cubic 
ipolation  (s  = 4),  regardless  of  the  fact  that  G-T's  accuracy  in 
time  is  only  of  second  order.  This  is  reflected  by  the  results  of 
Table  2. 
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